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ABSTRACT 

Aims. We are presenting here a study of the cold dust in the close environs of the ring nebula Cum 31. We aim at deriving the physical 
properties of the molecular gas and dust associated with the nebula, and investigating its correlation with the star formation in the 
region, that was probably triggered by the expansion of the ionization front against its environment. 

Methods. We make use of 870 pm emission data obtained with the Large APEX Bolometer Camera (LABOCA) to map the dust 
emission. The 870 pm emission provides an excellent probe of mass and density of dense molecular clouds. The obtained LABOCA 
image was compared to archival infrared, radio continuum, and optical images. 

Results. The 870 pm emission follows the 8 pm (Spitzer), 250 pm, and 500 pm (Herschel) emission distributions showing the classical 
morphology of a two dimensional projection of a spherical shell. We use the 870 pm and 250 pm images to identify 60 dust clumps in 
the collected layers of molecular gas using the Gausscluinps algorithm. The clumps have effective deconvolved radii between 0.16 
pc and 1.35 pc, masses between 70 Mq and 2800 Mq , and volume densities between 1.1 x 10^ cm“^ and ~ 2.04 x 10^ cm“^. The total 
mass of the clumps is ~ 37600 Mq . The dust temperature of the clumps is in the range from 21 K to 32 K, while inside the Hit region 
reaches 40 K. The clump mass distribution for the sample is well-fitted by a power law dN/(ffog(M/Mo) oc with a = 0.93 ± 
0.28. The slope differs from those obtained for the stellar IME in the solar neighborhood, suggesting that the clumps are not direct 
progenitors of single stars/protostars. The mass-radius relationship for the 41 clumps detected in the 870 pm emission shows that 
only 37% of them lie in or above the high-mass star formation threshold, most of them having candidate YSOs projected inside their 
limits. A comparison of the dynamical age of the Hit region with the fragmentation time, allowed us to conclude that the collect and 
collapse mechanism may be important for the star formation at the edge of Cum 31, although other processes may also be acting. The 
position of the identified young stellar objects in the region is also a strong indicator that the collect and collapse process is acting. 

Key words. ISM: molecules, Infrared: ISM, ISM: Hii regions, ISM:individual object: Cum 31, stars: star formation. 


1. Introduction 

There is now consensus that the formation of stars can be trig¬ 
gered by the action of Hit regions over its parental molecu¬ 
lar environment. Processes that sweep up and compress the 
gas and dust, like the co llect-and-collapse mechanism (C&C; 


for early stages of star formation (e.g. iRomero & Cappa 

2009 

Cappa et al. 

20091: Vasauez et aLll2012l: iDeharvens et al. 

2012 


tion process, especially the C&C mechanism, has been stud¬ 
ied extensively on the edges of many bubble-shaped Hit re- 


217, Sh2-90 (iDeharveng et al. 

120031 20081 

2009t IZavasno et aLl 

120061120101; iBrand et al.l 2011 

; ISamal et al. 

I 2 OI 4 I). as well as in 

seyeral suryeys le.g Thompson et al.ll2012; 

Simpson et al.ll2012^ 


Gum 31 is one of the many southern Galactic fJii regions 
that appear as a ring at infrared (IR) wavelengths (see the up¬ 
per panel of Fig. [3]l. The nebula, located at a distance of ~ 2.5 


kpc (lYonekura et al.ll200-5l). is considered to be part of the giant 


Carina Nebula Complex (iPreibisch et al.l 1^121) . and is ionized 
by the stars HD 92206A, B, an d C, which are members of the 
young stellar cl u ster NGC 3324 (I.Teffers et al.lll963 ^ Fort^[l97d 
WalbornI 119821: [Baumgardt et al.l I 2 OOOI: iMafz-Apellaniz et all 


lElmegreen & Lad^ 19771) and the radia tive driven implosion 
process (RDI: iLefloch & Lazarefil Il994t) . may favor the trig¬ 
gering of star formation. Dense molecular condensations (or 
clumps) lying at the border of Galactic Hit regions are then 


20041) . The simple morphology of Gum 31, along with the strong 


evidence of star formation in its environs (see below), make this 
object an excellent laboratory to investigate possible scenarios of 
triggered star formation. Cappa et al. (2008; hereafter CNAV08) 
analyzed the ionized, neutral, and molecular gas in the environs 
of Gum 31. Adopting a distance of 3 kpc, CNAV08 found an Hi 
shell of ~ 1500 Mq around the nebula with a systemic velocity 
of -23 km s ' expanding at ~ 11 km s *. The authors also found 
an associated molecular envelope (~ 1.1 x 10^ Mq ) at velocities 
between -27.2 km s“' and -14.0 km s“* that has been probably 
accumulated by the expansion of the ionized gas. The distribu¬ 
tion of the molecular and ionized gas, along with the emission of 
the 8 /im-MSX band A, which include significant emission from 
polyciclic aromatic hydrocarbons (PAHs), suggests that a pho¬ 
todissociation region (PDR) was created at the interface between 
the ionization front and the molecular cloud. Using MSX, IRAS, 
and 2MASS photometric data, CNAV08 also found a number 
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of infrared young stellar object (YSO) candidates around the 
nebula which suggests that star formation is active in the high 
density gas. As part of the Ce nsus of High - and M edium-mass 
Protostars (CHaMP) program. iBarnes et al.l (1201 ll) surveyed the 
southern sky to detect massive molecular clumps using Mopra 
(HPBW = 40") HCO^(l-O) observations. The authors identified 
19 molecular clumps in the region of Gum 31 (see Fig. 34 of that 
work). The systemic velocity of these clumps is in agreement 
with the velocity of the molecular gas reported by CNAV08. 
Ohlendorf et al. (2013; hereafter OPG13) used high sensitivity 
and spatial resolution WISE, Herschel, and Spitzer data to char¬ 
acterize the young stellar population surrounding the Hii region. 
The authors estimated a total population of ~ 5000 YSO can¬ 
didates in the region, many of them observed in compact clus¬ 
ters and located at the inner edge of the bubble. They concluded 
that probably the C&C and RDI processes are taking place si¬ 
multaneously at the border of the Hii region. The authors based 
their conclusion on the size of Gum 31, that agrees w ith what 
is expected for a C&C scenario (IWhitworth et al.iri994h . and the 
presence of YSOs in its rims, which give support to the C&C 
process, while the existence of Y SOs in the head of pill ars points 
to a RDI process. More recently. I Vazzano et al.l d20I4ll analyzed 
the distribution of the molecular gas and cold dust in the dense 
clump linked to IRAS 10361-5830, located in the southwestern 
edge of Gum 31. The authors found a dense molecular shell at ~ 
-25 km s * * centered at the IRAS position, where star formation 
is active. They also found a number of cold dust clumps associ¬ 
ated with the molecular shell. 

Although a study of the interstellar medium (ISM) and the 
young stellar populatio n around the whole Gum 31 nebula was 
performed by CNAV08,|B arnes et 10 (1201 Ih . and OPG13, ahigh 
spatial resolution analysis of the densest molecular gas adjacent 
to the edge of the Hii region still remains to be done. Low-J tran¬ 
sitions of '^CO and '^CO are commonly used to trace molecular 
clouds; however, they are often optically thick, thus probing only 
the outer layers of the cloud. Morover, also optically thin tran¬ 
sitions from C'^O and C'^O may fail to probe the inner dens¬ 
est molecular gas becaus e they freeze out onto dust grains and 
deple te at high densities dHernandez et al.ll201 it iGiannetti et al.l 
|2014|) . In this context, optically thin dust continuum emission 
in the (sub)millimeter range is one of the most reliable tracers 
of the dense molecular gas from which the stars form. It pro¬ 
vides a powerful tool to probe some basic physical properties of 
dense clouds (column density of molecular hydrogen, mass of 
the cloud, etc.) which are needed to unveil the physical condi¬ 
tions in regions where stars can form. 

When an Hii region expands, molecular gas and dust ac¬ 
cumulates behind the ionization front, forming shells of dense 
molecular gas surroun ding the ionized gas. With time these 
shells become massive dHosokawa & Inutsukall200^ and could 
contain cold dust that radiates in the (sub-)millimiter range. In 
the present study we analyze the distribution and physical prop¬ 
erties of the dense gas and dust adjacent to the ionization front 
of Gum 31 using 870 pm continuum emission obtained with 
LABOCA at APEX 12m telescope. The 870 pm image is com¬ 
pared to IR, radio continuum, and optical archival images in or¬ 
der to perform a complete multiwavelength study of the dust and 
gas associated with the nebula and to derive their physical prop¬ 
erties and conditions. Optically thin dust emission at 870 pm 
is usually dominated by the thermal emission from cold dust, 
which is contained in dense material (e.g. dense star-forming 
cores or filaments). 


2. Observations and complementary data 

The 870 pm continuum observations were carried out on 
October 2011 with the Large Apex BOlometer CAmera 
(LABOCAfl LABOCA is a 295-pixel bolometer array de- 
veloped by the Max-Planck-Institut ftir Radioastronomie 
(ISiringQ et al.ll200^ . The central frequency of the instrument is 
345 GHz and the bandwidth is 60 GHz. The beam size (HPBW) 
at 345 GHz is 1972 (or 0.23 pc at 2.5 kpc). Observations were 
made using the on-the-fly (OTP) mode to map an inverted 
“L”-shaped region of the nebula where the 8 pm emission is 
more intense. During the observations the amount of precip- 
itable water vapour (PWV) was between 0.16 mm and 0.24 
mm. Absolute flux calibrations were achieved through obser¬ 
vations of the planet Mars as primary calibrator for LABOCA, 
and the stars N207 HR and VY CMa as secondary calibrators. 
The uncertainty due to flux calibration was estimated to be ~ 
20%. The telescope focus and pointing were checked using the 
star rj Carinae. Observations were smoothed down to a beam of 
20'.'0 to obtain a final average rms noise of 50 mJy beam”'. The 
data were reduced using the Comprehensive Reduction U tiliPf 
for SHARC-2 software package (CRUSH-2fl (lKovacsl2008h fol¬ 
lowing the standard procedure. The 870 pm emission image of 
the Gum 31 is shown in Pig[I] 

The 870 pm data were complemented with several archival 
data sets; 

- Infrared images from the the Hersch el Infrared GALacti c 
(Hi-GAL) plane survey key program dMolinari et al.llTOloh . 
We used the 10 pm and 160//m maps produced in the High- 
Level 2.5, and 250pm and 500pm maps produced in the 
High-Level 3 (combined maps from Parallel Mode observa¬ 
tions done in both nominal and orthogonal scan directions) 
available in the Herschel Legacy as a standalone product. 
The Level 2.5 maps were produced with c orrectedmadmap 
in th e MADMap software application (ICantalupo et al.l 
1201 Ol) with point source artifacts corrected. The pixel sizes 
of the two PACS maps are 372 (for the 7 0 pm map) and 
475 (f or the 160 pm map), as suggested in iTraficante et al.l 
(1201 ih . The angular resolution for 70pm and 160pm bands 
are 10"and 1375, respectively. The Level 3 were produced 
from the Planck zero-point calibrated maps and the angular 
resolution for 250pm and 500pm maps are 18"and 36". 
The Herschel Interactive Processing Environment (HIPE 
vlfl Ott & Herschel Science Ground Segment Consortium! 
1201011 was used to handle the maps. 

- Infrared Spitzer images at 8.0 pm with a spatial resolution of 
~ 2", retrieved from the Galactic Leg acy Infrared Mid-Plan e 
Survey Extraordinaire (GLIMPSEn (iBeniamin et all2003h . 


’ APEX is a collaboration between the Max-Planck-Institut fiir 
Radioastronomie, the European Southern Observatory, and the Onsala 
Space Observatory 

^ http://www.submm.caltech.edu/~sharc/crush/index.html 
^ Herschel is an ESA space observatory with science in¬ 
struments provided by European-led Principal Investigator 
consortia and with important participation from NASA 
(http://www.cosmos.esa.int/web/herschel/science-archive I 

* HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS and SPIRE consortia members, see 
http://herschel.esac.esa.int/HerschelPeople.shtml 
^ http://sha.ipac.caltech.edu/applications/Spitzer/SHA 
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Fig. 1. 870 fim emission of Gum 31. The contours go from 150 mJy beam ' (~ 3cr) to 650 mJy beam ' in steps of 100 mJy beam 
and from 650 mJy beam^' in steps of 200 mJy beam The maximum flux density is 1860 mJy beam''. 


- Radiocontinuum data from the survey of t he Sydney 
University Molonglo Sky Survey (SUMSSS (iBock et al.l 
1 19991) at 843 MHz. 

- Optical data from the 2nd Digitized Sky Survey (red plate 
dMcLean et al.lbOOOl) . 


3. Observational results and discussion 

3.1. 870 fim emission and comparison with other 
waveiengths 

In Fig. [T] we show the 870 pm emission of the Gum 31 neb¬ 
ula. Despite the difference in spatial resolution, the 870 pm 
emission coarsely follows the CO(l-O) emission distribution de¬ 
picted in CNAV08 (see their Fig. 5). The emission at this wave- 
len gth shows the typical IR morphology of Galactic ring nebu¬ 
lae ^Ctorc^eFet2]|2Q^|20^|Deteryeng^et^i]20T3 with a 
very intense emission towards the border of the nebula and a very 

® http://www.astrop.physics.usyd.edu.au/cgi-bin/postage.pl 
’ http://skyview.gsfc.nasa.gov/cgi-bin/query.pl 


faint emission towards its center. From the figure, we note that 
the dust emission in the border is not uniformly distributed, and 
several condensations can be distinguished in the whole struc¬ 
ture. These condensations will be identified and analyzed in the 
next sections. 


The ring morphology in the IR emission of many Hii re¬ 
gions is a common feature i n the Galaxy (e.g. IWatson et all 
120081 : iDeharveng et alJ l2009i l2010t) and has originated some 
debate on whether these structures are flat (few par¬ 
secs) two-dimensional ring-like objects ( Beaumont^^Mlhamd 
20101) or three-dimensiona l bubbles dPeharveng et al. 120101 : 


Anderson et alj|201^ 12015]) . Although it is difficult to distin¬ 


guish between the two different scenarios, the observational 
properties suggest that the ring morphology of Gum 31 in the 
IR may be due to a projection effect, namely: 1) A weak struc¬ 
ture is observed in the 870 pm emission towards the center of 
the nebula along RA « lO'" 37™ 40^ (see Fig. [T]). This structure 
probably corresponds to IR emission of cold dust in front of the 
nebula, since its outer border shows an excellent spatial correla¬ 
tion with a region of optical absorption and 70 pm emission (see 
upper panel of Fig. [3]). The optically-absorbed region lies at RA 
> lO'' 37™ 45^, which lead us to conclude that probably there is 
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Observer 


Fig. 2. Sketch of Gum 31, which could explain in a qualitative way the IR, optical and radio continuum emission distribution. 


more 870 fim emission, arising from the front face of the bubble 
in the not-covered area with LABOCA. 2) The velocity interval 

of the h igh-density molecu lar gas adjacent to the nebula (-16 

km s~*: lBarnes et al.ll201lh . which indicates the presence of ex¬ 
panding motions that are not suitable to flat face-on structures. 3) 
The optical and radio continuum emission distribution inside the 
nebula (see below), which resembles a filled sphere of ionized 
gas inside the IR emission. 

In a qualitative manner, the morphology of Gum 31 (strong 
IR emission towards the border of the nebula and weak emission 
at the center) depicts the typical emission distribution of opti¬ 
cally thin emission of a bubble, scaling with the line-of-sight 
path length towards the edge of the shell. This effect is known as 
“limb brightening” (see Fig. |3. The weak 870 yum emission at 
the center may be due to gas/dust from the front face of the bub¬ 
ble with too low column density to be firmly detected. The limb 
brightening effect h as been used to explain the shape of man y 
IR ring bubbles (e.g. lWatson et al.l2008tlDeharveng et al.ll2010h . 
and could be also used to explain the spatial distribution of can¬ 
didate YSOs observed in the nebula, which seems to be more 
crowded at the edge of the bubble (see Fig. 8 of OPG13). High 
resolution molecular observations of the whole nebula would be 
useful to confirm the bubble morphology proposed above for 
Gum 31. 

In the upper panel of Fig. [3] we show the 870 yum emis¬ 
sion (white contours) and the SUMSS 843 MHz radio contin¬ 
uum emission (cyan contours) overlaid onto the optical (DSSR2) 
image (in green) and infrared (IRAC-GLIMPSE) 8 yum emis¬ 
sion (in red). At 8 yum, most of the emission originates in 
strong features of PAH molecules, which are considered to 
be good tra cers of warm UV-irradiated photodissociation re¬ 
gions (PDR; iHollenbach & Tielenslll99'^ . Since these complex 
mole cules are destroyed i nside the ionized gas of an Hit re¬ 

gion (ICesarskv et al.ll996l:IPovich et al.l2007l:lLebouteiller et all 
l2007h . they indicate the limits of the ionization front and de¬ 

lineate the boundaries of the bubble nebula. The 870 yum emis¬ 
sion follows the 8 yum emission distribution, which is bright in 
the western and southern regions of Gum 31. The 843 MHz 
radio continuum and optical emissions appear mostly confined 
within the IR bubble. Three emission peaks can be noticed in 
the radio continuum emission at RA, Dec.(J2000) = 10*37"'24*, 
-58°39'57", RA, Dec.(J2000) = 10*37'’T8", -58°34'50", and 


RA, Dec.(J2000) 10'’36”'59", -58°35'20". The last one is 

elongated and coincident with strong 870 yum emission, which 
strongly suggests that dense cold dust and molecular gas have 
been compressed and collected behind the ionization front due 
to the expansion of the Hii region. The other two radio contin¬ 
uum maxima are coincident with strong optical emission. 

The presence of intense emission in optical, 8 yum, radio con¬ 
tinuum, and 870 yum emission, to the west of the powering stars 
speaks in favor of the existence of a PDR at the interface be¬ 
tween the ionized and molecular gas, as suggested by CNAV08. 
On the other hand, the strong radio continuum, optical, and IR 
emission, observed in the western region of the bubble in com¬ 
parison with the emission at the eastern and northeastern edges 
suggests that the Hit region is expanding anisotropically, prob¬ 
ably due to the lack of dense molecular gas at the eastern and 
northeastern edges of the bubble. 

As can be seen in Fig. [2 a considerable deficit in the 8 
yum and 870 yum emission is visible around RA, Dec.(J2000) » 
10*36'"50^, -58°43'20", coincident with a lack of CO emission 
(CNAV08) and with the presence of some weak Ha and radio 
continuum emission located outside the border of the PDR. In 
Sect. 3.3 (Fig. |6li we will show that this region is correlated 
with a significant dust temperature enhancement. This could in¬ 
dicate that UV photons from the nebula could be leaking into 
the ISM through a discontinuity in the PDR. Two other interest¬ 
ing features can be noticed when comparing the 8 yum and 870 
yum emissions at RA, Dec.(J2000) 10''37'"03", -58°24'23"and 

RA, Dec.(J2000) ^ 10''37'”42", -58°27'00". The 870 yum struc¬ 
tures at these positions (which will be identified in Sect. 3.2 as 
clumps 8, 18, 19, and 35) appear projected onto deep spots of 
absorption at 8 pm emission, a typical feature of the so-called 
infrared dark clouds (IRDCs: [Egan et al.lll998l : [Rathborne et al.l 
I2006h . 

In the lower panel of Eig. [2 we show the 870 pm emission 
(green contours) overlaid on the 3 color image of the Herschel 70 
pm (in blue), 250 pm (in green), and 500 pm (in red) emissions. 
The 250 pm and 500 pm emission follow the 870 pm emission, 
although they seem to be a bit more extended. This means that 
the Herschel emission at high wavelengths is still good tracer 
of cold dust in Gum 31. A comparison between the 870 pm 
emission and the 250 pm and 500 pm emissions in the south¬ 
ern part of the nebula (Dec-58°46'00") shows that the area 
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Fig. 3. Upper panel: Composite image of Gum 31 and its environs. The IRAC-GLIMPSE emission at 8.0 pm is shown in red, while 
green indicates the DSSR2 optical emission (red plate). White contours show the 870 pm emission while cyan contours shows the 
radio continuum emission at 843 MHz. Black crosses at the center show the position of stars HD 92206A, B, and C. Lower panel: 
Composite image of the 870 pm dust continuum emission (green contours) overlaid onto the Herschel 70 pm (blue), 250 pm (green), 
and 500 pm (red) emissions. The contours corresponding to the 870 pm emission in both panels are the same as in FiglT] 


mapped with LABOCA almost covers the whole cold dust in the 
region. Since the emission at 250 pm and 500 pm is extremely 
faint at the eastern and north-eastern regi ons of Gum 31, where 
no CO(l-O) emission was detected either dYonekura et al.ll^05l 


CNAV08), we infer that there is little dense molecular gas/dust 
in the not-mapped area of the Hit region (see Sect. 3.3.4). This 
may indicate that Gum 31 has evolved in an inhomogeneous 
medium, with more material and higher densities towards its 
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western edge. The 70 /tm emission, which traces warm dust, is 
very intense along a very thin lane at the inner border of the 
PDR. This i s a common feature o bserved in many Galactic IR 
bubbles (e.g lAnderson et aDl2012l) . The 70 yum emission is also 
coincident with strong 870 yum emission in the southern region 
of the bubble, where the emission at 8 yum, 250 yum, and 500 
yum is also intense; this gives an idea about the complex nebula 
morphology. 

3.2. Identification of dust clumps at 870 pm 

As mentioned in Sect. 3.1, Figure[T]shows that the dust emission 
is not uniformly distributed, and several condensations can be 
distinguished in the whole structure. These condensations must 
be identified and their physical properties must be estimated if 
a study of the physical properties of the cold dust and molecu¬ 
lar gas in the whole bubble nebula is to be carried out. Different 
terminologies have been used in the last years to describe the 
substructure of a molecular cl oud. A hierarchical categorization 
was proposed by iBlitzi (Il993h . for structures of different scales 
in the interstellar medium, as: clouds, clumps, and cores. Clouds 
have masses >10 "^ Mq , and cores (regions where single stars or 
binary systems might be formed) have masses < 10 M© and sizes 
<0.1 pc. Clumps are then structures with intermediate masses 
and sizes from which stellar clusters can be formed. A simi- 
lar classification was also proposed by [Bergin & Tafallal (l2007h 
and IWilliams et al.l (l2000h designating individual overdensities 
inside a molecular cloud, which can be identified by different 
automated algorithms, as clumps. Given the characteristics men¬ 
tioned above and the physical properties derived for the conden¬ 
sations distinguished in the 870 pm emission, we will hereafter 
refer to these structures as “clumps”. 

We have identified the dust clumps in Gum 31 in a sys- 
tematic way with the commonly used Gaussc lumps algorithm 
(IStutzki & Guestenlll9^ iKrarner et al.lll998h . which is a task 
included in the GILDA^ package. Although Gaussclumps was 
originally written to work on 3-dimensional data cube, it can 
be modified to be applied to continuum images without modifi¬ 
cation of the code by adding two 2-dimensional empty planes 
to the origina l continuum imag e mimicking a 3-dimensional 
datacube (see iMotte et alJ l2003h . Gaussclumps uses a least- 
squares fitting procedure to decompose the image/cube into a se¬ 
ries of Gaussian-shaped clumps by iteratively substracting fitted 
clumps. The algorithm identifies the position, peak flux density, 
full width at half maximum (F WHM), and po s ition a ngle of the 
fitted sources. As suggested by iKramer et akl (Il998h . the “stiff¬ 
ness” parameters that control the fitting were all set to 1. The 
peak flux density threshold was set to 4cr. To avoid false detec¬ 
tions, we have trimmed the noisy borders of the image and added 
a wedge of null values to force Gaussclumps to identify clumps 
inside the image. To identify genuine clumps we have imposed 
that the deconvolved size of the clumps should be larger than 50 
% of the HPBW size (i.e. Rd > HPBW/2, or 0.12 pc at a distance 
of 2.5 kpc). As a result, 41 clumps were identified in the 870 pm 
emission image. 

In Table [T] we present all the clumps identified with 
Gaussclumps. The identification number of the clumps and 
their coordinates are indicated in Cols. 1 and 2, respectively. 
The spatial location, numerical identification, size, and orien¬ 
tation derived for each clump are depicted in Fig. |4] As can be 
seen from this figure, most of clumps satisfy a convincing vi¬ 
sual identification. In Col. 11 of Table [T] we indicate the HCO^ 


http://www.iram.fr/IRAMFR/GILDAS 
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Fig. 4. Clumps identified with Gaussclumps in the 870 pm 
emission image. 



Right Ascension (J2000) 

Fig. 5. Clumps identified with Gaussclumps algorithm in the 
250 pm emission image, over the not-covered region in 870 pm 
emission. 
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clumps reported bv iBarnes et al.l (1201 ih that are projected over 
the corre sponding 870 /rm clum p. Clumps 1, 6, and 3 were re¬ 
ported bv IVazzano et alJ (l2014l) . labeled in that work as Dl, D2, 
and D3, respectively (see Fig. 7 of that work). 

In Sect. 3.1, we pointed out that clumps 8, 18, 19, and 35 ap¬ 
pear projected onto deep spots of absorption at 8 /rm emission, 
which suggests that these clumps could be foreground dense ob¬ 
jects. These clumps, very bright at 250 jjm and 500 jjm but no 
emitters at 70 fim, appear quite isolated and projected slightly 
farther form the nebula, casting some doubts on their association 
with the rest of the clumps and their interaction with the ioniza¬ 
tion front. S ince their velocitie s are in agreement with the rest of 
the clumps (iBarnes et al.ll201 iT) we infer that these clumps could 
be a remnant of the parental cloud that remains undisturbed by 
the expansion of the nebula. 

With the aim of studying the north-eastern region of the Gum 
31 nebula, that was not observed with LABOCA (see Sect. 3.1), 
we make use of Herschel images at 250 //m since the emission 
at this wavelength is still a good tracer of cold dust. The iden- 
tihed clumps are presented in the lower part of Table [T] The 
number identihcation of the clumps are continued from those 
identihed at 870 yum. The spatial location, numerical identihca¬ 
tion, size, and orientation derived for each clump are depicted in 
Fig .|5] We have rejected from the sample two clumps around the 
position RA, Dec.(J2000) ^ 10'’38'"35", -58°42'30", since they 
likely are the 250 jum counterparts of clumps 9, 24, and 28. We 
also note that clumps 49,57, and 60 are projected over the center 
of the nebula and show a good spatial correlation with regions 
of optical absorption (see Fig.O which could indicate that these 
clumps are placed in front of the nebula. This clumps have low 
column and volume densities and could be related with the faint 
IR emision at the center of the nebula, giving more support to 
the bubble morphology proposed earlier for Gum 31. 

We are cautious about the clump identihcation described 
above, as well as the physical properties (see next section), since 
the IR nebula was proposed to be the result of a projection of 
a three dimensional bubble onto a two dimensional plane. This 
problem is, however, not only concerning to the case of Gum 31, 
but should be also addressed in many other molecular studies of 
Galactic bubble nebulae, where different molecular clumps (or 
condensations) were studied. 

3.3. Physical properties of the clumps 
3.3.1. Temperature 

In many works, it is assumed that (sub)m illimeter clumps/clouds 
would have a common temperature (e.g. lMookeriea et^l2004t 
iKirk et aDl2006t fLonez et al.ll2Ql ih . This is a logical assump¬ 
tion in the case of isolated clouds, athough it is not valid for 
regions like Gum 31, where different parts of the nebula are af¬ 
fected by different levels of radiation. Then, we have used the 
10 pm and 160jum Herschel images to obtain the dust tempera¬ 
ture of the clumps and of the whole Gum 31. To construct the 
Hux-ratio map, the 10 pm map was smoothed down to the an¬ 
gular resolution of the 160yum map. We divided the intensities 
of the map s by the recommend ed color-correction factors de¬ 
scribed bv iBernard et al.l (l20ir)h : 1.05 for the 10pm and 1.29 
for the 160/zm, values derived from a comparison with Planck 
data. The zero-level of the PACS archive data is unknown, which 
may cause problems to derive the background emission in each 
map. Here we adopted the background levels of 0.0023 Jy/px 
for 10 pm and 0.005 Jy/px for 160//m found bv iPreibi^h et alJ 
(l2012h . Then, the color-temperature map was constructed as the 


inverse function of the ratio map of Herschel 70 and 160jum 
color-and-background-correctedmaps, i.e., (see details 

in iPreibisch et alJl2012h . 

Assuming a dust emissivity following a power law /Cy oc v®, 
being j0 the spectral index of the thermal dust emission, in the op¬ 
tically thin thermal dust emission regime f(j) has the parametric 
form: 

_ 570 B{10,T) 110 f 

5160 B(160,7’) (l6oj 

where B{10, T) and B(160, T) are the Black Body Planck func¬ 
tion for a temperature T at the frequencies 70 pm and 160 pm, 
respectively. The pixel-to-pixel temperature was calc ulated as¬ 
suming a typical value [3-2 (see lPreibisch et afl2012l and refer¬ 
ences therein) on the whole Gum 31 region. The uncertainty in 
derived dust tempe ratures using this meth od was estimated to be 
about ~ 10 - 15 % (iPreibisch et al.ll20T^ . 

In Fig. |6] we show the color-temperature map obtained us¬ 
ing the method explained above. T he map shows a good cor- 
respondence wi t h that obtained in IPreibisch et al.l (1201 2h and 
lOhlendorf et al.l (1201 3l) . indicating temperatures < 28 K for 
clumps lying at the outer border of the bubble, while clumps 
projected over the inner region of the bubble (clumps 23, 28, 
and 60) show temperatures > 28 K, with the exception of clump 
49. Inside the bubble, close to the position of the stars of NGC 
3324 (where no 870 yum emission is detected), temperatures are 
above 35 K. The lowest derived dust temperature is 21 K, corre¬ 
sponding to clump 8. The most remarkable feature in Fig.|^is the 
excellent morphological correspondence of the 870 pm emission 
with darker regions in the temperature map, which confirms that 
we are dealing with the coldest and the densest molecular gas 
and dust that were collected behind the ionization front of the 
Hii region bubble. 

Temperatures of the clumps (Tjust, Col. 10 in TablellJ were 
simply estimated positioning the center of the htted Gaussian on 
the color temperature map. 


3.3.2. Mass 


The total (H 2 - 1 - dust) mass of the clumps (Mtot, Col. 7 in Table 
[TJ was calculated from their integrated 870 pm emission (^ 870 , 
Col. 4 in Table [T]i, as suming that the em ission is optically thin, 
using the equation of I Hildebrandl (1 198 3h : 


44(tot) — B 


I 870 I 


'‘"870 B870(7’dust) 


( 2 ) 


where R is the gas-to-dust ratio, d is the distance (adopted as 
2.5 kpc), /C 870 is the dust opacity per unit mass at 870 pm, as¬ 
sumed to be 1.0 cm^ g ' (e stimated for dust grains with thin 
ice mantles in cold clumps: lOssenkonf & Henning|[l994ll . and 
B 87 o(Tdust) is the Planck function for a temperature Tdust- For the 
case of the clumps detected in the Herschel 250 pm emission, 
we used a dust opacity a- 25 o = 12.1 cm^ g“', which is consistent 
with A-srn = 1.0 cm^ g~^ (assuming qc v^). Adopting R = 186 
(I.Tenkins[l2004l: IPraine et al.]l2007l : feeuther et^lioi ih . derived 
masses are in the range from 70 Mq to 2800 Mq and the total 
mass is ~ 37600 M© . This is a lower limit for the total mass in 
the nebula, since substantial 870 pm emission lying outside the 
clumps’s boundaries was not taken into account (see Fig. HI). In 
addition, some of the 870 pm emission towards the centre of the 
nebula (see Sect. 3.1) was probably not detected. We estimate a 
total mass uncertainty of about 50 -60 % arising mostly from a 
conservative distance uncertainty of ~ 20 % and a dust tempera¬ 
ture uncertainty of 15 %. 
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Table 1. Properties of identified clumps around the bubble Hit region Gum 31. Clumps 1 to 41 were identified at 870 jum while 
clumps 42 to 60 at 250 yum 


clump 

number 

RA, Decj 2 ooo 

m s o / //^ 

^peak 

(J^\ 

beam '' 

5v 

(Jy) 

^maj ^ ^min 
( X -) 

(pc) 

-^tot 

(Mo) 

n 

(cm“^) 

Ah2 

(lO^^cm^^) 

Tdust 

(K) 

HCO+ 
(BYF'ti ) 

1 

10 37 59, -58 46 34 

1.76 

23.41 

86.8x61.1 

0.42 

2050 

93500 

10.4 

25 

77b 

2 

10 37 51,-58 47 18 

1.71 

8.21 

45.7 X 41.8 

0.23 

770 

203800 

10.9 

24 

77a 

3 

10 38 07, -58 45 26 

1.59 

30.26 

113.9 X 66.5 

0.51 

2800 

72300 

9.9 

24 

77c 

4 

10 37 01,-58 33 59 

0.96 

16.62 

120.3 X 57.5 

0.48 

1450 

43800 

5.6 

25 

70a 

5 

10 37 38, -58 46 52 

0.92 

13.51 

87.5 X 66.7 

0.44 

1130 

43900 

5.2 

26 

77d 

6 

10 37 54, -58 45 57 

0.93 

4.65 

53.5 X 37.2 

0.23 

410 

103400 

5.5 

25 

77b 

7 

10 36 55, -58 36 34 

0.83 

6.11 

79.4 X 36.8 

0.29 

530 

70700 

4.9 

25 

70b 

8 

10 37 02, -58 24 27 

0.87 

7.21 

74.6 X 44.1 

0.32 

820 

84400 

6.7 

21 

67 

9 

10 38 37,-58 44 32 

0.79 

28.03 

130.1 X 108.7 

0.71 

2430 

23300 

4.6 

25 


10 

10 38 49, -58 47 10 

0.76 

32.98 

194.8 X 88.3 

0.78 

2640 

18900 

4.1 

26 


11 

10 38 24, -58 45 36 

0.73 

18.59 

132.8 X 76.1 

0.59 

1640 

26900 

4.4 

25 


12 

10 36 54, -58 37 52 

0.56 

4.52 

56.7 X 56.1 

0.31 

390 

40600 

3.2 

25 

70b 

13 

10 36 51,-58 28 47 

0.56 

3.65 

69.3 X 37.1 

0.27 

370 

60800 

3.9 

22 

66 

14 

10 36 56, -58 31 30 

0.52 

6.10 

96.8 X 48.2 

0.39 

580 

33600 

3.3 

23 

69 

15 

10 38 53, -58 44 58 

0.54 

3.56 

105.0 X 24.8 

0.23 

290 

76200 

3.0 

26 


16 

10 37 58, -58 44 01 

0.51 

5.84 

97.4 X 46.2 

0.38 

480 

29600 

2.9 

26 


17 

10 38 36, -58 49 20 

0.51 

6.95 

127.8 X 42.3 

0.41 

540 

26100 

2.7 

27 


18 

10 37 43, -58 27 49 

0.48 

4.93 

77.8 X 52.7 

0.36 

470 

32500 

3.1 

23 

72 

19 

10 37 45, -58 26 22 

0.48 

3.75 

58.2 x 53.6 

0.31 

370 

40600 

3.2 

23 

72 

20 

10 38 34, -58 47 48 

0.48 

16.91 

130.1 X 106.5 

0.70 

1290 

12800 

2.5 

28 


21 

10 37 25, -58 46 48 

0.46 

8.71 

97.5 X 77.4 

0.51 

700 

17900 

2.5 

27 


22 

10 38 52, -58 49 03 

0.50 

1.94 

57.2 X 26.8 

0.18 

150 

77100 

2.6 

27 


23 

10 37 46, -58 37 42 

0.40 

2.30 

69.2 X 33.2 

0.25 

140 

30100 

1.7 

32 


24 

10 38 50, -58 43 33 

0.38 

2.68 

71.5 X 39.1 

0.29 

200 

27900 

1.9 

28 


25 

10 37 35, -58 42 59 

0.36 

4.86 

91.1 X 58.9 

0.42 

390 

17400 

1.9 

27 


26 

10 37 27, -58 48 17 

0.34 

5.48 

99.0 X 63.7 

0.46 

470 

16300 

2.0 

25 


27 

10 37 46, -58 43 43 

0.34 

3.44 

78.1 X 51.8 

0.36 

270 

19100 

1.8 

27 


28 

10 38 46, -58 41 44 

0.33 

1.05 

46.5 X 26.8 

0.16 

70 

52700 

1.5 

28 


29 

10 36 43, -58 38 07 

0.33 

4.79 

90.5 X 63.0 

0.44 

390 

15600 

1.8 

26 


30 

10 37 46, -58 47 53 

0.34 

2.37 

61.3 X 44.5 

0.29 

170 

23500 

1.6 

29 


31 

10 38 09, -58 46 37 

0.34 

1.47 

44.9 X 38.0 

0.21 

120 

38600 

1.8 

27 


32 

10 37 02,-58 32 31 

0.32 

2.68 

98.3 X 33.7 

0.30 

240 

27900 

1.9 

24 


33 

10 38 08, -58 44 20 

0.35 

2.35 

83.4 X 32.2 

0.27 

180 

30900 

1.8 

27 


34 

10 38 49, -58 44 47 

0.33 

2.09 

69.1 X 36.1 

0.27 

170 

30200 

1.9 

26 


35 

10 37 37, -58 26 50 

0.31 

1.74 

59.0 X 37.4 

0.25 

160 

34700 

2.0 

24 

72 

36 

10 36 56, -58 30 07 

0.31 

2.69 

64.8 X 53.1 

0.33 

240 

22100 

1.9 

25 


37 

10 37 46, -58 45 53 

0.33 

3.22 

64.1 X 60.4 

0.35 

240 

17700 

1.6 

28 


38 

10 36 49, -58 47 51 

0.29 

6.04 

117.7 X 70.0 

0.53 

520 

11700 

1.7 

25 


39 

10 36 58, -58 35 11 

0.30 

2.75 

78.0 X 46.2 

0.34 

220 

19200 

1.6 

28 


40 

10 37 39, -58 30 35 

0.28 

12.56 

187.0 X 94.5 

0.79 

990 

6800 

1.5 

27 


41 

10 37 09, -58 33 07 

0.29 

1.73 

72.8 X 32.0 

0.25 

130 

27600 

1.5 

28 



42 

10 39 05, -58 42 40 

8.94 

329.38 

147.3 X 80.9 

0.65 

550 

6800 

1.2 

25 

43 

10 39 28, -58 42 36 

7.31 

464.92 

162.3 X 126.9 

0.86 

810 

4400 

1.1 

24 

44 

10 39 10, -58 40 25 

6.91 

270.21 

127.7 X 99.2 

0.67 

470 

5300 

1.0 

24 

45 

10 38 34, -58 33 21 

6.21 

153.44 

101.8 X 78.5 

0.53 

240 

5600 

0.8 

25 

46 

10 38 47, -58 35 20 

6.04 

170.64 

103.5 X 88.4 

0.56 

270 

5100 

0.8 

25 

47 

10 38 11, -58 30 46 

5.69 

128.66 

90.4 X 81.0 

0.50 

200 

5300 

0.7 

25 

48 

10 38 36, -58 27 15 

5.94 

394.71 

156.4 X 137.6 

0.88 

680 

3400 

0.8 

24 

49 

10 38 07, -58 39 45 

5.45 

129.48 

94.8 X 81.1 

0.52 

180 

4400 

0.6 

27 

50 

10 38 18, -58 28 34 

5.59 

397.16 

208.4 X 110.5 

0.91 

520 

2400 

0.6 

27 

51 

10 39 15, -58 29 52 

5.44 

846.49 

242.6 X 207.4 

1.35 

1390 

1900 

0.7 

25 

52 

10 39 06, -58 38 37 

5.28 

149.84 

120.4 X 76.3 

0.56 

230 

4300 

0.7 

26 

53 

10 38 40, -58 31 16 

5.29 

368.83 

178.7 X 126.2 

0.90 

510 

2400 

0.6 

27 

54 

10 39 45, -58 41 16 

5.33 

292.23 

147.5 X 120.4 

0.79 

520 

3450 

0.8 

24 

55 

10 39 21,-58 25 54 

5.20 

312.22 

160.3 X 121.3 

0.83 

600 

3500 

0.8 

23 

56 

10 38 42, -58 24 47 

5.09 

285.67 

184.5 X 184.5 

0.80 

610 

4000 

0.9 

23 

57 

10 38 17, -58 41 04 

4.79 

290.35 

176.9 X 111.0 

0.84 

360 

2100 

0.5 

28 

58 

10 38 49, -58 33 44 

4.69 

175.54 

131.5 X 92.0 

0.65 

270 

3300 

0.6 

25 

59 

10 39 22, -58 36 27 

4.23 

597.59 

251.0 X 182.3 

1.29 

890 

1400 

0.5 

26 

60 

10 38 27, -58 37 46 

3.75 

542.53 

263.4 X 177.8 

1.30 

670 

1100 

0.4 

28 


Notes. From Table 4 and Figure 34 of iBarnes et al.l ll201lh . 
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Fig. 6. Dust temperature map (in color scale) derived from Herschel emission at 70 and 160 //m. The color-temperature scale is 
on the right. The 870 fim emission is shown in green contours and the 250 jim emission is shown in white contours. The 870 fim 
contours are the same as in Fig.[T]and the 250 jjm contours go from 410 mJy beam * (~ 4cr) to 590 mJy beam ' in steps of 50 mJy 
beam”', and from 700 mJy beam”' in steps of 100 mJy beam”*. 


3.3.3. Size 

The effective deconvolved radius of the clumps (/?d, Col. 7 in 
Table [T]i was calculated as 

Rd = X (3) 

where 0maj and ^min are the major and minor deconvolved 
FWHM of the clump (Col. 5 in Table [TJ. These values were ob¬ 
tained as: 

^maj ^ tlx “ ^HPBW and ^jnin ^ dy “ ^HPBW ’ 

being dx and dy the major and minor FWHM sizes derived by 
the Gaussclumps algorithm, and 6hpbw the final resolution of 
the image (20" for 870 yum and 18" for 250 jum). Clump radii 
are found to lie between ll'.'l and 6572 (0.16 pc and 0.79 pc, 
respectively, at a distance of 2.5 kpc) for 870 yum, while for 250 
yum clumps radii are between 4173 to 11174 (0.5 pc to 1.35 pc 
at a distance of 2.5 kpc). The difference could reside in that the 


emission at 250 yum may be tracing warmer and less dense dust 
than the 870 yum emission making its emission more extended 
(see Fig. |3]l. This could also impact in the determination of the 
volume densities of the clumps, which are slightly lower than 
those determined from the 870 yum emission (see next section). 

3.3.4. Densities 

The average volume density of each clump in. Col. 8 in Table[T]i 
was derived, assuming a spherical geometry, as: 

n - -r- (5) 

4/3 7T niii 

where fi is the mean molecular weight (assumed to be equal to 
2.8 after allowance of a relative helium abundance of 25% by 
mass), and mu is the mass of the hydrogen atom. Volume den¬ 
sities of the clumps are between 1.1 x 10^ cm”^ and 2.03 x 10^ 
cm”^. The average density in the sample is 2.9 x 10^ cm”^. 
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The beam-averaged column density of the clumps , Col. 
9 in Table [T]l was calculated using 


= R 


^peak 


f^beam OTh fi87o(Tdust) 


(6) 


(lHildebrandlll983h . where fpeak is the 870 jjm emission peak in¬ 
tensity (Col. 3 in Table [1]) and fibeam is the beam solid angle 
(;7r ^hpbw / lii(2)). As we did when determining the mass, we 

assumed R = 186. Considering a lower limit for the solid an¬ 
gle corresponding to the HPBW size, an intensity peak limit of 
3cr (150 mJy beam^^ for 870 fim, and 400 mJy beam^' for 250 
fj.m), and a maximum dust temperature for the clumps of 32 K 
(see Table[T]|, the minimum beam-averaged column densities de¬ 
tectable at 870 fim and 250 /rm are 0.63 x 10^^ cm“^ and 0.32 x 
10^^ cm“^, respectively. 

Column densities derived for all the identified clumps are 
between 0.38 x 10^^ cm“^ and 10.85 x 10^^ cm“^. The average 
column density of the sample is ~ 2.4 x 10^^ cm“^, almost twice 
the value obtained by CNAV08 for the molecular shell associated 
with Gum 31. This is not surprising given that we are probing the 
molecular gas intimately associated with the ionization front of 
the Hit region, which is more affected by the expansion of the 
nebula. We keep in mind that this difference could also be due to 
different angular resolutions and abundance variations. We also 
note that volume densities of the 870 fim clumps are, in average, 
larger than those of 250 fim clumps by a factor of ~ 10. The same 
occurs with column densities, which are larger by a factor of ~ 
5. A possible explanation for this difference is an evolution of 
the Hit region in an inhomogeneous cloud. The faint 8 fim, 250 
fim, and 500 fim emission towards the eastern and northeastern 
regions of Gum 31 supports this scenario (see Sect. 3.1). On the 
other hand, as suggested in Sect. 3.3.3, the 250 fim emission 
could be tracing warmer and less dense dust than the 870 fim 
emission. A combination of both factors is also possible. 

The physical properties derived for the dust clumps in 
Gum 31 (radius, density and mass) are similar to those 
found in IR continuu m clumps detected in other Galactic star¬ 
forming regions (e. g. iMookeriea et al.ll2004 iLonez et'aDl201 it 
iBeuther et al.ll201 ih . 


3.4. Clump mass distribution 

Using the mass obtained for each clump, we derived the clump 
mass distribution (CMD) in the Gum 31 region. The CMD is 
plotted as £/N/£/log(M/Mo) versus mass. The first term is ap¬ 
proximated as the number of clumps in each mass interval (N) 
and the second is approximated by the logarithmic mass inter¬ 
val, Alog(M/Mo). The error was determined considering only 
the statistical Poisson uncertainty, VAN. We also fitted the his¬ 
tograms with a power law of the form dN/dlogiM/Mo) oc 
where a is the spectral mass index. A CMD can also be plotted as 
dNIdiMfM^), which has a power-law dependence £/N/£/(M/M©) 
oc where T = cr H- 1 . 

To overcome fitting artifacts due to different bin sizes 
(Alog(M/Mo)), we derived spectral indexes for different bin 
sizes, ranging from 0.25 to 0.45 in steps of 0.01. For each bin 
logarithmic bin size we fitted the corresponding spectral mass 
index. In all the cases, a turnover in the distribution can be seen 
at bins between 300 M© and 500 M©, which is likely caused by 
the low-mass incompleteness. To make sure that the fit is not af¬ 
fected by incompleteness effects at lower masses, we use a lower 
cut-off of 400 M© for the fit. In Fig. |7]we show an example of a 
CMD derived for Alog(M/M©) = 0.31. The procedure resulted 


in 21 fits, obtaining spectral indexes a in the range 0.44 - 1.35, 
with a final weighted mean spectral index a = 0.93 + 0.28. 

Although making a comparison of the spectral mass indexes 
is a complex task due to the many assumptions authors made, 
especially when extracting clumps and deriving their masses, 
the average spectral mass index derived for the dust clumps 
in Gum 31 is of the same order of others derived for molec- 


ular clouds from dust continuum (e.g 

. Mookeriea et all 12004 

ILonez et al.l201 ItiBeuther et al.l201 It 

VIiettinenl2012t) and from 

molecular line CO observations (e.g. 

31itzl 1993;lKrame 

et al.l 

1998; Simon et al.ll2001h iHever et al.l 

200 It Simon et al.l 

20061; 


Wong et alJl20(m 



Fig. 7. Example of a CMD for Alog(M/M©) = 0.31. Red trian¬ 
gles represent the central mass of each bin and error bars rep¬ 
resent the standard deviation of a Poisson distribution ( VAN). 
Blue lines represent the best fit above the 400 M© limit (dotted 
vertical line). 


The slope of a CMD is u sually compare d to that of the initial 
stellar mass function (IMF. Salneter 1955 ) as a way to asses a 
direct relation clump-star (iMotte et al.l 19981: iLada et akl 120071 : 
iLonez et"^ 1201 ih . Equal slopes would indicate that clumps 
form stars directly. The weighted average slope derived for Gum 
31 is shallower than Salpeter’s value (a - 1.35 or F = 2.35) 
and other indexe s between 2.1 and 2.5 in the uppe r mass regime 
(above ~ 1 M© ; lKrour)all200ll : I Andre et'^l2010l) implying that 
the clumps are not direct progenitors of individual stars. This 
means that other processes, such as fragmentation (mainly of 
the most massive clumps) are necessary to determine the initial 
stellar masses. We keep in mind, however, that observations with 
higher spatial resolution, resolving core structures o f ~ 0.1 pc in 
size, could result in a spectral index of ~2.35 (e. g. iMotte et al.l 
Il998h . We also emphasize that the radii of the clumps are be¬ 
tween 0.16 pc to as large as 1.35 pc in size (most of them hav¬ 
ing sizes above 0.4 pc), which are about the size of star-forming 
regions, rather than individual star-forming cores (typically be¬ 
tween 0.01 pc and 0.1 pc in size). Probably, smaller clumps (or 
cores), may not be correctly separated with the resolution of the 
LABOCA data, although this is just speculation. Since these hy¬ 
pothetical smaller cores can not be identified in Herschel im¬ 
ages, higher spatial resolution images (HPBW < 8", equivalent 
to sizes < 0.1 pc at a distance of 2.5 kpc) of the whole nebula 
would be useful to shed some light on this issue. 
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Fig. 8. Candidate YSOs projected over the 870 yum emission. WISE Class I and Class II candidates (OPG13) are marked by 
red and green diamonds. MSX massive young stellar object (MYSO) and compact Hit region candidates (CNAV08) are marked 
by magenta and yellow circles, respectively, while 2MASS candidate YSOs (CNAV08) are indicated by blue crosses. Molecular 
clumps identibed in Sect. 3 are indicated by blue ellipses. 



3.5. Star formation evidence 


According to lLada & Ladal (l2003h and iMotte et al.l (l2003h the ra¬ 
dius and mass required to form stellar clusters are about 0.5 pc 
to 1 pc and 100 Mq to 1000 Mq , respectively, which are similar 
to those obtained with Gaussclumps in Gum 31. The clumps 
identibed with Gaussclumps in the 870 yum emission are then 
good candidates to form clusters of stars, rather than individ¬ 
ual stars. As pointed out in Sect.l, an exhaustive search of YSO 
candidates was carried out by OPG13 using high sensitivity and 
high spatial resolution IR data. Although some contamination 
from the fore- and background sources is present, OPG13 con¬ 
sidered the WISE-selected sample of YSO candidates to be the 
more reliable, bnding 661 candidates (207 Class I and 454 Class 
II sources) of solar to intermediate mass. Since most of the can¬ 
didate YSOs detected by OPG13 and CNAV08 appear projected 
onto the area mapped with LABOCA, and the gas and dust in the 
eastern and northeastern regions of Gum 31 is much less dense 


than the southern and western regions, we will concenttate our 
analysis only on the 870 yum emission. 

In Eig.[8]we show the location of the WISE YSO candidates 
reported by OPG13 that appear projected over the 870 yum emis¬ 
sion. The location of the YSOs were cautiously determined us¬ 
ing a composite image of 8.0 /urn, 4.5 yum, and 3.6 yum emissions 
as a reference. We have also included the MSX and 2MASS 
candidate YSOs identibed by CNAV08. Erom a visual inspec¬ 
tion of Eig. |8] we can roughly distinguish bve clusters of YSO 
candidates projected onto the submillimeter clumps. Two small 
compact clusters are approximately centered at RA, Dec.(J2000) 
^ 10 A 37 moo^, -58°34'47" and RA, Dec.(J2000) ^ 10''36'"54^ 
-58°37'26". The brst cluster is composed by 10 protostellar 
members and is projected onto clumps 4, 32, and 39. The second 
cluster is composed by 13 members and is projected over clumps 
7, 12, and 29. In the southern region of the bubble, three clus¬ 
ters can be discerned. The brst one, approximately centered at 
RA, Dec.(J2000) 10'G8'"42", -58°45'45", contains approxi¬ 

mately 31 members projected over clumps 9, 10, 11, 15, 20, and 
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Fig. 9. Mass-radius relationship. Filled red traingles indicate 870 fim clumps with candidate YSOs projected inside the clumps’s 
limits determined by Gaussclumps, while green open triangles indicate 870 yum starless clumps. Magenta open squares indi- 
cate the 250 um clumps. T he blue solid line shows the scaled mass-radius th reshold for high mass star formation determined by 
iKauffmann & Pillail ( |2010|) . The rescaled mass-radius relationship derived by lUrauhart et al.l (1201 3h is shown by the dotted black 
line. Upper and lower solid red lines show surface densities of 1 g cm“^ and 0.024 g cm“^, respectively. Dotted green lines indicate 
volume densities at 10^ cm“^, 10^ cm“^, 10^ cm“^ , and 10® cm“^. 


34. The second one, located at RA, Dec.(J2000) 10'’38'"02", 

-58°45'42", consists of ~ 29 members projected over clumps 1, 
2, 3, 6, 16, 30, and 33, while the third one is approximately cen¬ 
tered at RA, Dec.(J2000) 10'’37'"13", -58°47'10", and is com¬ 

posed approximately by 18 members projected only over clump 
21 . 

In Fig. 0 we show the mass-radius relationship for the 41 
clumps identihed in the 870 //m emission, discriminating be¬ 
tween clumps that have candidate YSOs projected inside their 
limits and clumps without signposts of star formation. For com¬ 
parison purposes only, we have also included the 19 clumps 
ide ntified in the 250 yum emis sion. _ 

IKauffmann et al.l (l2010allbll and IKauffmann & Pill^ (1201 Ol) 
investigated the mass radius relationship of nearby molecular 
clouds and found the empirical relationship M [Mq] < 1972 x 
(^eff [pc])' as the threshold to avoid high mass (M ~ 10 Mq ) 
star formation (hereafter the K&P threshold). Note that when 
deriving the ir relationship, Kauffm a nn et al. reduced the dust 
opacities of lOssenkopf & Hennind (Il994ll by a factor of 1.5. 
Since this correction has not been applied here, w e have rescaled 
the re lationship given in the original work (cf. lUrauhart et alJ 
1201 3ll . The rescaling was also calculated considering the opaci¬ 
ties and gas-to-dust ratio used in this work (/cg 7 o = 1.0 cm^ g“* 
and R = 186; see Sect. 3.3). The K&P threshold is shown in 
Fig. |9] For the sake of comparison, we also show the empir¬ 
ical mass-size relationship obtained by lUrauhart et alJ (1201 3l) 


for intermediate- and high-mass star formation, also rescaled 
to /cgTO = 10 cm^ g * and R = 186. The upper and lower red 
solid lines in Fig. |9] depict the constant surface densities of 1 g 
cm“^ and 0.024 g cm“^, respectively. The surface density of 1 
gcnG^ (47 87 Mp/pc^) wa s advo cated by iKrumholz & McKeel 
(120081) and iMcKee & Taiil (1200^ as criterion to avoid exces¬ 
sive fragmentation and to allow massive star formation. The sur¬ 
face density of 0.024 g c m~^ represents the average surface den- 
sity threshold de rived by iLada et al.l (1201 Oh (116 Mopc^^) and 
iHeiderman et al.l (1201 Oh (129 Mq pc ^), for efficient star forma¬ 
tion. 

As can be seen from Fig. |3 all clumps lie above the lower 
surface density limit of 0.024 g cm“^. Of the 41 clumps identi¬ 
hed in the 870 /rm emission, only 15 (37 %) lie above the K&P 
limit, with 12 clumps (clumps 1, 2, 3,4,6, 7, 8, 9, 10,15, 11, and 
20) having candidate YSOs projected inside their limits. Only 
clumps 5, 13, and 14 are starless clumps. As mentioned i n Sect . 
3.2, clumps 1 and 3 were also found bv IVazzano et al.l (l2014h 
(labeled in that work as D1 and D3, respectively) and were asso¬ 
ciated by the authors to a dense molecular shell linked to several 
infrared sources. Clump 1 is coincident with the MSX sources 
0286.3747-00.2630 and 0286.3773-00.2563, which were clas- 
sihed as a compact Hii region (CHii) and massive young stellar 
object (MYSO) candidates, respectively (sources 20 and 21 in 
CNAV08). Two WISE Class I and Class II candidates identihed 
by OP013 are also projected inside the boundaries of this clump 
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(see Fig. |8]l. This clump have an HCO^ counterpart (BYF 77b; 
IBarnes et al.ll201 ih which confirms the existence of high density 
molecular gas. Regarding clump 3, five WISE candidate YSOs 
identified by OPG13 are projected onto it. This clump also have 
an HCO^ counterpart (BYF 77c). Clump 2 also have an HCO^ 
counterpart (BYF 77a) and is coincident with the MSX source 
G286.3579-00.2933 (identified by CNAV08 as a MYSO candi¬ 
date), the 2MASS candidate YSO source 10375219-5847133, 
and two WISE candidate YSOs. For the case of clump 7, a very 
crowded spot of candidate YSOs appear projected inside its lim¬ 
its. The peak intensity of this clump is almost coincident with 
the position of the MSX CHii candidate G286.1626-00.1877 and 
the 2MASS sources 10365396-5836293 and 10365749-5836366 
(sources 27 and 29 in CNAV08). Three WISE candidate YSOs 
(OPG13) also appear in side the clump. Thi s clump has an HCO^ 
counterpart (BYF 70b: lBarnes et al.ll201 ih . Clumps 6 and 8 are 
very dense and also have HCO^ counterparts (BYF 77b and BYF 
67, respectively), although only have two and one WISE candi¬ 
date YSOs projected inside their limits. Regarding clump 4, it 
has an HCO^ counterpart (BYF 70a) and appears projected onto 
a number of WISE sources although many of them seems to be 
located at the inner region of the bubble, outside the molecu¬ 
lar gas (see Fig. 8 of OPG13). Clumps 9, 10, and 11 are less 
dense (they do not have HCO^ counterparts), but a considerably 
number of WISE candidate YSOs are seen projected inside their 
limits. 

Several clumps with candidate YSOs projected on them lie 
below the K&P mass-radius relation limit, namely: clumps 12, 
16, 20, 21, 29, 36, and 39. Probably, the most interesting clumps 
of this sample are 12 and 39, since they are located in a region 
where the star formation activity seems to be very intense (see 
previous paragraph). Clump 12 has three WISE candidate YSOs 
projected inside, while clump 39 has six. Both, clumps 12 and 39 
have HCO^ counterparts (BYF 70b and BYF 70a, respectively). 

3.6. Triggered star formation scenario 

As mentioned before, the morphological characteristics of Gum 
31 make this nebula a perfect object to investigate possible sce¬ 
narios of star formation. To confirm whether the fragmentation 
of the collected layer of molecular gas and dust have occurred 
(as sugge sted by PPG 13 and in S ect. 3.4), we use the analytical 
model of I Whitworth et all (Il994h for the case of an expanding 
Hii region which predicts the time at w hich the fragmentatio n 
occurs. We have followed the analysis of IZavagno et al.l (l2006h . 
deriving the time when the fragmentation may have occurred 
(Ifiag), and the size of the Hii region at ffrag (Ruag), given by 

ffrag = 1.56flo.2'/"n3^^"e(H%-'/" lO'^ yr (7) 




Fig. 10. Upper panel: Fragmentation time (red curves) and dy¬ 
namical time (green curve), for R - 5.5 pc as a function of no = 
0.2,0.4, and 0.6 km s“'. Lower panel: radius of fragmentation as 
a function of no for the same values than in the upper panel (red 
curves). The green curve represents the radius of the Hii region 
(5.5 pc). 

where R is the observed radius of the nebula, ahii is the speed of 
sound in the ioniz ed medium, and R p is the radius of the initial 
Stromgren sphere (I Strom grenll 193 9h given by 


Arag = 5.8 no.2^/" n 3 '’/" e(H% 




pc 


where ao .2 is the isothermal sound speed in the compressed layer, 
in units of 0.2 km s“' (aJO.! km s“'), n^ is surrounding homo¬ 
geneous infinite medium into which the Hii region expands, in 
units of 1000 cm“^ (uq/IOOO cm“^), and ^(H”)^^ is the number 
of ionizing photons, in units of 10"^® s“' ((2(H°)/10"^® s“'). In or¬ 
der to determine f^ag, in Fig. [10] (upper panel) we compare Eq.|2l 
with the stan dard evolutionary model fo r expanding Hii regions 
developed bv iDvson & WilliamsI (Il997h . 


R = Ro 


, Gyn 

1 -H OhII X 

Ro 


4/7 


(9) 


Ro[«o, e(H°)] = 3.15 X 10-'’ nf/^ e(H°)'/3 pc (10) 

For the case of NGC3324, two 06.5V and one 08.5 are the 
main sources of UV photons (see Sect.l), which provide a total 
ionizing flux of ~ 1.7 x 10''® s ' (ISmith & Brooksll2007h . Since 
about 50 % of the U V photons co uld be absorbed by interstellar 
dust in an Hii region (Ilnouell20'()i1) . we adopt a final value (2(H°) 
~0.85x 10''® s-'. 

The upper panel of Fig.fTOlshows ffrag and Gyn as function of 
the ini tial ambient density hq- As pointed out in lWhitworth et all 
(11994 . a, = 0.2 km s ' is likely a lower limit for the sound 
of speed in the collected layer, since both turbulence and ex¬ 
tra heating from intense sub-Fyman-continuum photons leaking 
from the Hii region could increase this value. Then, in the fol- 
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lowing we will adopt a range Us - 0.2 - 0.6 km s ' for the col¬ 
lected layer. For the ionized medium we adopt a value ahii - 10 
km s *. Since fdyn is required to be higher than fuag for fragmen¬ 
tation, the initial ambient density in the region where Gum 31 
was formed must have been in a range no ~ 950 - 1800 cm“^. 
A slightly larger range (no ~ 1100 - 2500 cm“^) is derived by 
comparing R with the radius of fragmentation (Fig. [10] lower 
panel). We have compared these values with the initial ambi¬ 
ent density derived using the mass estimates from our 870 yum 
data and mass estimates from CNAV08. If we assume that all 
the material observed in the nebula next to the ionization front 
was swept up during its expansion, an estimate of no for Gum 31 
can be obtained by assuming that all the gas (either ionized or 
neutral) was initia lly located in a sphere of radius equal to that of 
the Hn region (e.gIZavagno et al.ll2007llCichowolski et al.l2009l: 

lAnderson et al.ll2015l) . For the observed radius of the bubble, R, 

we consider the angular radius of ~ 7(5 derived from the radio 
continuum emission at 843 MFIz (CNAV08) which yields to R 
= 5.5 pc at a distance of 2.5 kpc. Considering a total molecular 
mass of 35800 Mq , derived from the 870 yum and 250 yum emis¬ 
sion without taking into account the mass of clumps 8, 18, 19, 
and 35 (see Sect. 3.2) and an ionized gas mass of 2430 Mq (at 
d - 2.5 kpc) associated to the nebula (3500 Mq at d - 3 kpc; 
CNAV08), and R - 5.5 pc, we derive an initial ambient density 
of tiQ ~ 2200 cm“^. This is a lower limit for no since a con¬ 
siderable amount of 870 yum emission outside the clumps and 
from the centre of the nebula was not considered (see Sects. 3.2 
and 3.2.2). Furthermore, we do not know if the original parental 
cloud extended to the current radius of the nebula. The derived 
value of no is higher than the density limits required for the frag¬ 
mentation according to the Whitworth et al.’s scenario, which 
indicates that the C&C process might be important in the col¬ 
lected layers of molec ular gas. _ 

As pointed out by IZavagno et al.l (l2006h . the estimates pre¬ 
sented above are rough, since the Whitworth et al’s model as¬ 
sumes expansion into an homogeneous medium and many un¬ 
certainties remain concerning to the uniformity of the medium 
where Gum 31 has evolved (see Sect. 3.1). Also, we agree with 
OPG13 that the presence of cores, bright rims and pillars at the 
edge of the bubble suggest that other processes, like the inter¬ 
action of the io nization front with pre-ex isting molecular con¬ 
densation (RDI; iLefloch & Lazarefflll994h or small-scale Jeans 
grav itational instabilities in the swept-up layers of molecular gas 
(e.g. lPomares et aO2009l) could be also active in some regions. 

4. Summary 

Using LABOCA 870 yum continuum emission and IR and opti¬ 
cal archival data, we have carried out a study of the gas and dust 
adjacent to the Gum 31 nebula. The main results can be summa¬ 
rized as follows: 

1. The morphology observed in the 870 yum emission for 
Gum 31 is typical of many Galactic IR ring nebulae. The 
emission is strong towards the border of the Hit region and 
faint towards its center, depicting the typical emission of a 
two dimensional projection of a spherical bubble (i.e. limb 
brightening effect). The IR emission distribution coincides 
with that of the molecular gas and suggests that the cold dust 
and molecular material have been collected behind the ion¬ 
ization front due to the expansion of the nebula. The intense 
optical, 8 yum, radio continuum, and 870 yum emissions at the 
west of the powering stars suggest the existence of a PDR 
at the interface between the ionized and molecular gas. The 


high intensity in radio continuum, optical, and IR emissions, 
in the western region of the bubble, as well as the low 
intensities at the eastern region suggest that the Hii region 
is expanding anisotropically, probably due to the formation 
of the bubble in an inhomogeneous parental molecular cloud. 

2. We have used the Gaussclumps algorithm to identify 
dust clumps adjacent to the nebula. After applying some 
rejection criteria, we were left with 41 clumps in the 870 
yum image. To detect and investigate de dust clumps in the 
eastern and northestern region of the nebula, which was 
not covered in the 870 yum emission data, we made use 
of an Herschel image at 250 yum and we identified 19 clumps. 

3. Adopting a distance of 2.5 kpc, the total (gas -H dust) clump 
masses are in the range from 70 Mq to 2800 Mq and the 
total mass is about 37600 Mq . The deconvolved linear 
radius derived for the clumps are between 0.12 pc and 1.3 
pc. Volume and column densities are in the range from ~ 
1.1 X 10^ cm“^ to ~ 2.04 x 10^ cm“^ and 0.4 x 10^^ cm“^ 
to 10.4 X 10^^ cm“^, respectively. Dust clump temperatures 
obtained from high resolution Herschel data at 70 yum and 
160 yum are between 21 K and 32 K, while dust temperatures 
at the center of the bubble, close to the ionizing stars of 
NGC 3324 are as high as ~ 40 K. 

4. We have derived the clump mass distributions for Gum 31, 

approximated as AN/Alog(M/Mo), using different bin sizes, 
and we fitted a power law (7N/£/log(M/Mo) for each 

bin size. Obtained spectral indexes are in the range 0.44 
- 1.35, and the weighted mean index is a = 0.93 + 0.28. 
The derived spectral index is of the order of other indexes 
derived in many works fro m continuum and molecular line 
data for other regions (e.g. iLonez et"^ 1201 it IWong et ^ 
l2008h . Although the resolution of LABOCA could not be 
sufficient to make a comparison, the difference between the 
average slope of the CMDs with the slopes of the stellar 
IMF found in the literature suggests that the clumps are not 
direct progenitors of single stars. 

5. The mass-radius relationship for the 41 clumps identified 
in 870 yum shows that all the clumps lie above the lower 
surface density threshold for massive star formation of 0.024 
g cm“^, while only 15 clumps lie in or above the K&P limit. 
Of the sample of clumps lying above the K&P threshold, 
only three of them do not have candidate YSOs projected 
inside its their limits. 

6. After analyzing the dynamical age of the Hit region and 
the fragmentation time and radius of the molecular/dust 
shell surrounding the Hit region, we confirmed that the 
collect and collapse process may indeed be important 
in the collected layers of gas at the edge of the bubble. 
The presence of YSO candidates detected by OPG13 and 
CNAV08 projected onto the 870 yum emission adjacent to 
the ionization front is a strong indication in favor of this 
scenario. We keep in mind, however, that other process like 
radiative-driven implosion or small-scale Jeans gravitational 
instabilities in the swept-up layers of molecular gas, could 
also be active in specific regions. 
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